2D Radial Distribution Function of Silicene 
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Silicene is the counterpart of graphene and its potential applications as a part of the current elec- 
tronics, based in silicon, make it a very important system to study. We perform molecular dynamics 
simulations and analyze the structure of a two dimensional array of Si atoms by means of the radial 
distribution function, at different temperatures and densities. As a first approach, the 2D Lennard- 
Jones potential is used and two sets of parameters are tested. We find that the radial distribution 
function does not change with the parameters and resembles the corresponding to the (111) surface 
of the FCC structure. The liquid phase appears at very high temperatures, suggesting a very stable 
system in the solid phase. 
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El siliceno es la contraparte del grafeno y sus potenciales aplicaciones en la electronica actual, basada 
en el silicio, lo hacen un objeto de estudio muy import ante. Realizamos simulaciones de dinamica 
molecular de un sistema bidimensional formado por atomos de Si, a diferentes temperaturas y 
densidades, analizando la estructura por medio de la funcion de distribucion radial. Como primera 
aproximacion, usamos el potencial de Lennard- Jones en 2D. Se encontro que utilizando dos conjuntos 
de parametros diferentes, la funcion de distribucion radial no cambia y se asemeja a la del piano 
(111) de la estructura FCC. La fase liquida aparece a muy altas temperaturas, sugiriendo un sistema 
muy estable en la fase solida. 

Descriptores: Siliceno; Funcion de distribucion radial; Dinamica molecular 



PACS numbers: 61.46.-w; 61.48.-c;68.65.-k 

I. INTRODUCTION 

During the last two decades, honeycomb-structured 
materials became very important in nanoscience. Their 
unique orbital symmetry gives rise to exceptional prop- 
erties in these quasi- ID carbon-based systems, such as 
fullerenes, nanotubes, graphene and its nanoribbons [l]- 
[3]. Among them, graphene is doubtless the most stud- 
ied at present time, due to its potential applications 
and the fact that is the first bidimensional material of 
one-atom thickness [4-6 . On the other hand, the inter- 
est on fabrication, characterization and study of silicon 
nanostructures keeps continously expanding. The pos- 
sibility of having a graphene-like structure, the silicene, 
has been studied theoretically already for several years 
[THTO]. From an electronic point of view, silicene could 
be equivalent to graphene, having the advantage of an 
easy integration to the present electronics, which is based 
on bulk silicon. Moreover, its borders do not react with 
oxygen [11 . Nevertheless, it has a drawback. Less versa- 
tile than carbon, silicon hardly hibridizes to sp'^ bonds. 
Therefore, syntesis and growth of silicene is extremely 
hard [12 . Only recently, silicon structures resembling 
graphene, such as self- assembled silicene nanoribbons [6] 
and silicene sheets deposited on silver crystals [13] have 
been reported. Growth on a ZrB2 substrate has been suc- 
cesful [H]. For the hexagonal lattice, the lattice constant 
a for Si and Ge, as reported by Lebegue and Eriksson [12] 
are 3.860 A and 4.034 A respectively, which are larger 



than the one for graphene (2.46 A) since these atoms 
have larger radii. Cahangirov et al [7] showed that the 
bonding distance Si-Si (Ge-Ge) is 2.25 A(2.38 A). 

With four sp electrons in valence band. Si and Ge are 
chemically very similar to C, in spite of this they behave 
very differently. The difference in the chemistry exhib- 
ited by carbon and silicon can be traced to the difference 
in their tt bonding capabilities. First, the energy differ- 
ence between the valence s and p orbitals for carbon is 
about twice that for silicon (Si: Esp-Es^ = 5.66 eV, C: 
E2p-E2s=10.60 eV) [15 ); as a result, silicon tends to uti- 
lize all three of its valence p orbitals, resulting in sp^ 
hybridization while, in contrast the relatively large hy- 
bridization energy for carbon implies that this will "ac- 
tivate" one valence p orbital at a time, as requirements 
the bonding situation, giving rise, in turn, to sp, sp'^ , and 
sp^ hybridizations. Second, the tt-tt overlap decreases by 
roughly an order of magnitude in going from carbon to 
silicon due to the significant increase in atomic distance, 
resulting in much weaker tt bonding for silicon in com- 
parison with that for carbon. Hence, Si=Si are in gen- 
eral much weaker than C=C bonds. [16]. Fagan et al. 
|17 establish the theoretical similarities and differences 
between Si and C nanotubes. For the silicon nanotubes 
(SiNTs) studied, they obtained a cohesive energy value of 
0.83 eV/atom, which is higher than the total energy per 
atom for the diamond-like structure. Considering that 
the cohesive energy for the Si bulk in diamond structure 
is 4.63 eV/atom, the energies for the studied nanotubes 
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are only 82% of the bulk. Comparing with carbon nan- 
otubes (CNTs) that have around of 99% of the cohesive 
energy they would have in perfect crystalline graphite, 
these results help to understand the diffilculty to pro- 
duce the SiNTs and, especially, silicene. Zhang et al. 
[16] also studied silicon nanotubes and their results sug- 
gest that silicon nanotubes can in principle be formed, 
however, the energy-minimized SiNTs adopt a severely 
puckered structure (with a corrugated surface) with Si- 
Si distances ranging from 1.85 to 2.25 A. 

Some silicene properties have been predicted theoreti- 
cally to be very similar to their corresponding in graphene 
[H [13] . Ab-initio calculations revealed that silicene clus- 
ters can be used in Field-Effect Transistors and hydro- 
gene storage, and their electronic properties have been 
studied using molecular dynamics [18 . Ince and Erkoc 
(IE) determined the structure of silicene nanoribbons of 
different widths and lenghts [iF. 

In this paper we report results on structural proper- 
ties of a bidimensional array of Si atoms. By means of 
molecular dynamics, we determine the radial distribution 
function for different temperatures and densities. We an- 
alyze the effect of the two-body potential considering a 
Lennard- Jones potential with two sets of parameters, the 
IE parameters and those suggested by Stillinger and We- 
ber (SW) for bulk Si [20 . Our calculations represent a 
description of free-standing silicene, since ab-initio stud- 
ies have shown that two dimensional honeycombe struc- 
tures for Si and Ge are stable [7 . Our model provides 
insight to understanding the structural properties of sil- 
icene. Taking into account potentials developed for sili- 
con or interactions with the substrate for the case of ad- 
sorbed nanosheets will provide an improved description. 
In Section 2 we briefly describe our calculation method 
and the used parameters. In Section 3 the results are 
analized and finally. Conclusions close this report. 



II. METHODOLOGY 

We perform molecular 2D dynamics simulations. Si 
atoms interact via Lennard- Jones (LJ) potential 

The first term is a positive (repulsive) short-range in- 
teraction, related to the electrostatic repulsion. The sec- 
ond term, a negative contribution, is the Van der Waals 
potential. The parameter e (a) corresponds to the energy 
(spatial range) scale. LJ potential works well for rare 
gases and poorly for materials where many-body effects 
are important, like metals or semiconductors. Neverthe- 
less, it can give a rough description and has been used as 
the two-body part in Si nanostructures [191 EI] • Details 
of our approach can be found in [22l [23] . 

We compare calculations performed with two sets of 
parameters shown in Table |l] The first set takes <j as the 



distance in the Si- Si dimer [19] while SW determined the 
values for crystalline Si at K [20]. They are alike, with 
the IE well being deeper and wider. 



(7 IE 


^lE 


0.2295 nm 


2.817 eV 
(4.5134x10-^^ J) 


crsw 




0.20951 nm 


2.168 eV 
(3.4738x10-^^ J) 



TABLE L LJ potential parameters, IE [l9] and SW [20]. 

We focus on the radial distribution function (RDF), 
which describes how the atoms in a system radially ar- 
range around each other. The RDF gives information 
about the average structure of disordered molecular sys- 
tems such as liquids. It is also of fundamental importance 
in thermodynamics, because some macroscopic thermo- 
dynamic quantities can be calculated using the RDF, for 
example the pressure and the energy. The RDF is a func- 
tion of the radial distance r, defined as 




where V is the volume, N is the number of atoms, Vij is 
the position vector of atom j respect to the ith atom, and 
the brackets indicate average over all atoms. This func- 
tion g{r) has characteristic shapes for different phases 
[2311M]. 

III. RESULTS 

Simulations with 512 atoms at different coverages (den- 
sities) and temperatures were performed. Densities were 
taken in dimensionless units (pj^)= 0.459 (0.604) 

-diamond Si-, 0.592 (0.778) and 0.690 (0.907), the latter 
ones higher than liquid silicon density [25 . The cho- 
sen temperatures, TJ^ (Tj£;)= 0.0119 (0.0092), 0.0670 
(0.0515) and 0.1987 (0.1530), correspond to normal con- 
ditions, bulk melting point and a higher temperature, 
respectively, since the melting point of 2D structures is 
known to increase considerably [26 . Table [ll] shows the 
values of density and temperature in real units. LJ po- 
tential tends to generate a FCC structure, and its cleav- 
age plane is (111) [27 . Besides, LJ tends to generate 
this kind of lattice. Thus, structures related to FCC are 
expected. 

Figs. [1] and [2] show the RDF's for IE and SW param- 
eters. Here r* = r/cr. At p^j^ {p%w)= ^-^^^ (0.459) and 
TJ^; (TJ^)= 0.0092 (0.0119) peaks are narrow, a typ- 
ical characteristic of the crystalline structure for solids. 
These RDF's coincide with the RDF in 2D obtained by 
Rodriguez et al. for Argon, but with a slight shift to 
the right. We see that even at very high temperatures the 
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Density 


Temperature 


/9, gr/cm"^ 


Psw {p}e) 


T, K 




2.33 


0.459 (0.604) 


300 


0.0119 (0.0092) 


3.00 


0.592 (0.778) 


1685 


0.0670 (0.0515) 


3.50 


0.690 (0.907) 


5000 


0.1987 (0.1530) 



TABLE IT Density and temperature, in real units, for poten- 
tial parameters IE [19] and SW [20]. 



peaks slightly broaden while the positions of the peaks 
and the overall shape do not change, indicating a very 
high stability of the system. The RDF at T}^ (TJ^) = 
0.1530 (0.1987) resembles the radial distribution func- 
tion for atoms in layer 1 of the (111) surface obtained by 
Broughton and Gilmer [28 , when they studied the struc- 
tures of atomic layers in the crystal- vapor interfaces for 
FCC structures by molecular dynamics. This FCC sur- 
face is characteristic for 2D honeycomb structures, with 
the difference that it has an extra atom at the center 
of each hexagon. It is worth mentioning that compared 
to graphene, the FDR's obtained are similar, showing a 
difference in the position of the peaks [29, 30 , with a 
scaling which depends on the value of the bond length, 
and increased peak heights due to contributions of the 
central atoms. 

In the RDF at T}^ {T^sw)= 0.1530 (0.1987) the peak 
heights diminish, because the short-range interaction be- 
tween atoms is being lost and the system approaches to 
the liquid phase. Similar results are found for the inter- 
mediate value of density, not shown here. 

At P*iE {Psw)= 0-907 (0.690), we obtained the RDF's 
shown in Figs |3] and [4j For the SW case, the RDF have 
similar behavior as described above. We can conclude 
that there is a negligible change with increasing density. 

However, for the IE case, the last peaks are almost 
lost when the density increases; this loss occurs roughly 
around r*= 3 for T*= 0.0092 and about r*= 4 for 
T*=0.0515, being this behavior more noticiable at T*= 
0.0092. At T*= 0.1530, the RDF does not show a difer- 
ent behavior to that presented at p*= 0.604. 

Finally, Figs. [5]and[6]show RDF's at very high temper- 
atures, and clearly their shapes correspond to the liquid 
phase, indicating that melting ocurrs at high tempera- 
ture. These high temperatures may be due to the re- 
striction on atoms to remain in two-dimensional motion, 
so it takes more energy to break sp'^ bonds, which are 
much less common for Si 112 . 



IV. CONCLUSIONS 

We have performed 2D molecular dynamics simula- 
tions to study the radial distribution function of silicene. 
A Lennard- Jones potential was considered and two sets of 
parameters were used. The RDF's obtained correspond 
to a two-dimensional close-paked hexagonal lattice, i. e. 
the plane (111) of the FCC structure. No significant dif- 
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FIG. 1. Radial Distribution Function, IE parameters and 
0.604 which corresponds to crystalline (diamond) silicon 
density. 
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FIG. 2. Radial Distribution Function, SW parameters and 
p*= 0.459, corresponding to crystalline (diamond) silicon den- 
sity. 
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FIG. 3. Radial Distribution Function, IE parameters and p^— 
0.907, higher than the corresponding to hquid siUcon density. 
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FIG. 4. Radial Distribution Function, SW parameters and 
p*= 0.690, higher than the corresponding to liquid silicon 
density. 



ferences between results obtained with the two sets of 
parameters were found. Our results suggest a very high 
melting temperature for the silicene. More detailed stud- 
ies with more appropiate potentials for covalent semicon- 
ductors are in progress. 
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FIG. 5. High temperature radial distribution function calcu- 
lated with IE parameters. Density and temperatures as 
indicated. A typical shape for a liquid appears. 




FIG. 6. High temperature radial distribution function calcu- 
lated with SW parameters. Density p* and temperatures as 
indicated. A typical shape for a liquid appears. 
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